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We have used the notion of the constituent quark model of nucleon, where a 
constituent quark carries its own internal structure, and applied it to deter- 
mine nuclear structure functions ratios. It is found that the description of ex- 
perimental data require the inclusion of strong shadowing effect for x < 0.01. 
■ Using the idea of vector meson dominance model and other ingredients this 

CN ■ effect is calculated in the context of the constituent quark model. It is rather 

striking that the constituent quark model, used here, gives a good account of 

CD ' the data for a wide range of atomic mass number from A = 4 to A = 204. 
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I. INTRODUCTION 



The measurement of the nuclear structure function, F^(x,Q 2 ), in Deep Inelastic Scat- 
tering (DIS) indicates that the structure function of a bound nucleon is different from a 
free nucleon. Very precise data on the proton structure function from HERA |J [0 |J ||] 
and NMC || which cover a wide range of kinematics, both in x and Q 2 , sets a stringent 
demand on the understanding of nucleon structure function in the nuclear environment. 
The difference between the parton distributions in the free and the bound nucleon, that is, 
f±(x, Q 2 ) 7^ fi(x, Q 2 ) is attributed to the nuclear effects. These effects are characterized in 

A p 

the structure function ratio 

p ^2_0^ Q 2 ) m 

Fj>(x,Q 2 ) [i) 

There are a host of interpretations that describe the DIS data on the nuclear targets pio- 
neered by the European Muon Collaboration (EMC effect) and it is believed that various 
mechanisms are responsible for the behavior of R in various kinematical intervals: 
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i) The Shadowing effect at x < 0.1 where R < 1. 

ii) Ant i- shadowing effect at 0.1 < x < 0.3 where i? > 1. 

iii) The EMC effect for 0.3 < x < 0.7 where R < 1. 

iv) The Fermi motion for x > 0.7 where i? > 1. 

For a detailed study of these mechanisms and other theoretical models sea Arneode |J. 
Currently, there are fairly accurate data on the structure function ratio R = -jinrf^m that 
has been obtained by different experimental groups for a variety of nuclei ranging in mass 
number from A = 4 to A = 208. These data, now, permit us to investigate the x and Q 2 
dependence of the F^(x, Q 2 ). One feature of the data, almost universal, is the shape of the 
^-dependence of the _F 2 A and is basically the same for every nuclei, thus making it rather 
easy to fit in a wide range of models with totally different underlying assumptions. Our 
interest here is to find a description for the nucleon structure which can be extended to the 
nuclei; yet maintaining some universal features common to both. To this end, our starting 
point is the study of the internal structure of the constituent quark which is common to 
both nucleon and nuclei. Its internal structure is universal and governed by the annihilation 
and pair production in Quantum Chromodynamics. We assert that at high enough Q 2 , it is 
the structure of this constituent quark that is being probed, whereas, at low Q 2 it behaves 
as a structure-less valence quark. Such a model has been applied successfully to the proton 
structure for the entire range of kinematics in x and Q 2 |§ . Encouraged by this results, we 
think that the nuclear environment can only alter the distribution of the constituent quark, 
leaving its internal structure intact. The fact that valence quark distribution in nuclei re- 
ceives a distortion as opposed to its distribution in a free nucleon is very similar to the 
nucleon swelling and the relaxation of the confinement of the constituent quark within the 
boundaries of a nucleon embedded in a nucleus. To put it differently, a nucleon in a nucleus 
will relax its boundaries and will occupy broader space, preserving the universality of the 
constituent quark of the relaxed nucleon. None of the number of the constituent quarks, its 
internal structure, and the average momentum of the constituent will change. Thus, here we 
will be dealing with a collection of structure-full U, D-type constituent quarks participating 
in DIS. 



II. FORMALISM 

To proceed quantitatively, we need to employ a model which is compatible with the deep 
inelastic scattering data on a free nucleon and then we impose the restrictions due to nuclear 
effects for the bound nucleons.The picture that comes in mind is the so-called Valon model 
of R.C.Hwa 0. In this model the nucleon is considered as a system of bound constituent 
quarks, which themselves have structure. The bound state problem is a nonperturbative 
effect which is taken into account in the distribution of constituent quarks. The structure 
of constituent quark is produced perturbatively and is free of bound state problem. In such 
a picture the structure function of a free nucleon is the convolution of constituent quark 
distribution in a nucleon and the structure function of the constituent quark itself: 

F 2 N (x,Q 2 ) = 24 J dyGu(y)f c (x,Q 2 ) + e 2 D J dyGn(y)f c (x,Q 2 ) (2) 
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where Gu_(y) and Go(y) are the distribution of U and D-type constituents in a nucleon 

and ff(z, Q 2 ) is the structure function of the constituent calculated in QCD. At sufficiently 
high Q 2 , ff(z, Q 2 ) can be described accurately in the Leading-Order results in QCD. They 
are also calculated in the Next-to- Leading Order [8]; but here we will restrict ourselves to 
the leading order, because, we are mainly interested in the application of the model and it 
describes the nuclear structure functions sufficiently accurately. The moments of ff(z,Q 2 ) 
can be expressed in terms of the evolution parameter: 

s = In — 4_ (3) 

m AS 

where Qq = 0.233^|^ 2 and A = 0.345^p 2 are scale parameters determined from the data. 
The moments of the singlet and non-singlet constituent quark structure function in the 
leading order solution of the renormalization group equation are given as[9]: 

M NS (n,Q 2 )=exp(-d NS s), . . 

M s (n, Q 2 ) = |(1 + p) exp(-d + s) + \{l - p) exp(-cL s), { > 

The anomalous dimensions, d' s , and other associated parameters are : 

p = (d NS - d gg )/A, 
A = d + - d_ = [{d NS - dgg) 2 + 4:d g Qd Qg } 1/2 , 
d *s = [1 - j], 

j _ -2 2+n+n 2 

U 9Q 3nb n{n' 2 -iy (r\ 
j -J 2+n+rfi \°) 

a Q9 2nb n(n+l){n+2) ' 

U 99 Ttbl 12 ^ nin-l) ^ (n+l)(n+2) 18 ^2 jh 

d± = l[d NS + dgg ± A], 
b = (33 - 2/)/12tt, 

Since these moments, M(n,s), are known, using inverse Mellin transformation technique, 
we can obtain the distributions of various components in the constituent quark; namely, 
valence, sea quarks, and gluon distributions, P v ( sea <9) . They are as follows: 

P^(z,Q 2 ) = ^-J c dnz- n+1 M v{s) g(n,s), 

M v (n,s) = M NS (n,s), 
M s (n, s) = (2f)-\M s (n, s) - M NS (n, s)), 1 ' 

M g (n,s) = Mg Q (n,s), 

and M g Q(n,s) is the quark-to-gluon evolution function given by: 

M gQ (n, S) = A^dgQlexpt-d+s) - exp(-cLs)], (7) 

To account for the SU(2) asymmetry of the nucleon sea which is evident form the violation 
of the Gottfried sum rule, we follow the same procedure as in and replace the distribution 
of d-quark in the sea by: 
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Pd/c(z, Q 2 ) = -^Z^Pu/ciz, Q 2 ), (8) 

The calculation of the moments given in eq.(4) is simple. Instead of exhibiting the moment 
distribution, we present the results in parametric form. That is, for every s we fit the 
moments by a form give below: 

P»(z, Q 2 ) = a v (Q 2 ) zMQ 2 ) (1 _ z yAQ 2 ) y 

P6(z,Q 2 )=Ela s AQ 2 )(i-z) b ^ Q2 \ (9) 

P£(z,Q 2 )=Zla g SQ 2 )(l-z)K(Q 2 ) 

for the valence, sea quarks, and gluon respectively. The dependence of the coefficients dj, bj 
and Cj on Q 2 , or rather on s is given in the appendix. To complete the calculation of , we 
also need to specify constituent quark distributions, Gu_ and Gd. The simplest approach 
to evaluate these distributions is to write down the inclusive momentum distribution of the 
constituent quark in a free nucleon. We will limit ourselves to proton here, these inclusive 
distributions can be written as: 

GWd (2/1,2/2,2/3) =« 2/1 2/2 + 2/2 + 2/3-1) (10) 

where a = 0.65 and b = 0.35 are two free parameters and j/i, 2/2 refer to [/-type and 2/3 refers 
to the Z?-type constituents. By double integration over the unspecified variable we get the 
exclusive distributions: 



Gu(y) = [B(a + 1, a + b + 2)]" 1 y a (1 - y) a + h + l 
Gn(y) = [B(b + 1, 2a + 2)]- V (1 - V? a+l 



(11) 



where B(x,y) is the Euler Beta function. The above distributions satisfy also the following 
normalization conditions: 

f 1 dyGu(y)= [ 1 dyGo(y) = l (12) 
Jo p Jo p 

In figure 1, we present the result for Ff(x, Q 2 ) for several values of x and Q 2 , its agreement 
with the experimental data in a wide range of kinematics is rather good. 
Having determined the nucleon structure function, next we extend the same procedure to 
the nuclear medium. As we mentioned earlier, in the nuclear medium the constituent quark 
distribution has to be modified, but its structure function, or rather parton distribution in 
a constituent will not be affected by the nuclear environment. The modification to Gc(y) is 
achieved by introduction of a distortion factor, 5. Hence, we write the inclusive momentum 
distribution as follows : 

G'uuoiy) = a( yi y 2 r +s y b 3 +5 5( yi + y 2 + y 3 - 1) (13) 

which leads to : 

G',(2/)=a^ (l-^+2W 

G f o{y)=a' y b+& {l-y) 2a+2&+1 (14) 
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with 



a = [Beta(a + 5 + 1, a + b + 25 + 2)]" 1 
a' = [Beta(b + 5 + 1, 2a + 25 + 2)] 



(15) 



we note that 5 is a function of atomic number, A. In Ref. || a probabilistic argument is 
used to find 5. Here we do generalize the same result as follows : 

0.0104 - 0.014 In A + 0.014(lnA) 2 
~ 1 - 0.83 In A + 0.295(lnv4) 2 - 0.024(ln A) 3 ^ 

So far we have determined the constituent quark distribution in a bound nucleon. We need 
to find its distribution in a nucleus. For that matter, however, one needs also to know the 
nucleon distribution in the nucleus, we can use the Fermi distribution and arrive at the 
following result : 



d 3 k p{k) 8{z - 1 - (17) 

k F lVl n 



where k + = k° + k z and k° = are the light cone coordinates and p(k) = j^pr@{kF — \k\) is 
the nucleon density in the nucleus. So, straightforward integration would yield the nucleon 
distribution in a nucleus. However, this also means that the bound nucleon is nearly on 
shell and therefore, the nuclear binding effect is neglected. To be more realistic, and include 
this effect, we replace z — 1 by z — T)a where t\a = 1 — and Mt = — Ba is the 
effective mass of the nucleon. With these fine tunings we arrive at the following function for 
the nucleon distribution in a nucleus : 

= !f(ft - (* - ^) 2 ) ->VA-£(z{VA+£ (18) 

4>n(z) = ; otherwise 

The values of the Fermi momentum, kp, and the binding energy, Ba, are given by the 
requirements of nuclear physics. Going one step further, we carry the Fermi motion to the 
constituent quark level. This is done easily by renormalization of the distribution function 
give in Eq.(14) this leads to the following constituent quark distribution per proton in a 
nucleus : 

Gi(y) « {-)G'u{-) + ^ 2 (-^)( V-(-) U (19) 
p r] p rj 10 az z z p z 




where A = -fe- . Having collected all the ingredients, now we are in a position to evaluate the 

N 

nuclear structure function ratios, R = J D , . Notice that all of our discussions about the 
constituent quark distribution are pertinent to proton. Finally, taking the neutron excess in 
the nucleus into consideration, we get : 
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F 2 A (x,Q 2 ) = (±){F 2A (x,Q 2 ) - ~(iV- Z)[F?(x,Q 2 ) -F|(*,Q 2 )]} (21) 

where F 2 (x, Q 2 ) and F£{x, Q 2 ) are the free proton and free neutron structure functions, 
respectively, and F 2 a(x, Q 2 ) is obtained using Eq.(l). The results of the model is presented 
in figure 2 for a variety of nuclei ranging from A = 4 to A = 204. As it is apparent from 
Fig. (2), we see that the model calculation is reasonable only down to x > 10 _1 . For smaller 
values of x, x < 10 _1 an extra ingredient is required. An obvious and natural place to look 
for the dynamics of small x would be the shadowing effects, where the growth of parton 
densities are suppressed due to recombination and annihilation of the partons. There are 
some attempts in the literature to calculate the shadowing corrections, pioneered by Qiu 



12 1 . However we found that such a parameterization is inadequate in describing the data. 
Following Ref. [12] we modify the shadowing correction of Qiu in two ways: (1) we take the 
general idea of vector meson dominance which describes low Q 2 virtual photon interaction, 
to model the approximate scaling of shadowing by pomeron exchange. It turns out that the 
experimental data on F 2p (x, Q 2 ) up to x < 0.05 can be well described by a sub-asymptotic 



pomeron with the effective anomalous dimension going down logarithmically with x ||13||.(2) 



The magnitude of the shadowing effect in such models is sensitive to the value of vector 



meson-nucleon cross section. This sensitivity is explored by G. Show [14] and it is described 
by a quantity : 

Art = (22) 

A Aa 1 N 

This ration is particularly manifested for heavy nuclei and reflects the ratio of nuclear radius 
R to the mean free path, A, of the vector meson inside the nucleus; which also depends on 
the nuclear density . That is, the shadowing effect at small x values increases with increasing 
nuclear radius and density. Utilizing these modifications to the original work of Qiu [12], 
finally we present the following modification of parton distributions in a nucleus : 



Pc/a( z ' Q 2 ) - Rs(%i Q 2 , A eff )P£, /N (z, Q 2 ;i ^3) 



where superscript s(v) stands for sea (valence) components and R s (x,Q 2 , A e ff) is the shad- 
owing factor affecting parton distributions of the constituent quark. These distributions in 
the nuclear medium are compared to the corresponding distribution in a free nucleon. The 
parameterized form of this factor is : 

R s (x,Q 2 ,A eff ) = R s (Q 2 ,A eff )(0M + 0MQA° e f f )x°- 5 ln(x)) (24) 

And Aeff is given by: 

A eff = -10.97 - 0.704A + .0042A 2 - 0.00019A 25 + 8.63A - 5 (25) 

In equation (24) R S (Q 2 , A e ff) is evaluated by the n* h -moment equations for the parton 
distribution in the sea of a constituent quark [[HJ: 

< F S C/N {Q 2 ) >„=< F° /N {Ql) > n K- q (Q 2 )+ < F 9 C/N {QD > n ^ 9 (Q 2 )+ 
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where < F(Q 2 ) > n = /d z n ~ 2 P(z,Q 2 )dz. Obviously, nuclear shadowing do not affect the 
distribution of valence quarks in the constituent quarks. If the momentum loss owing to 
the recombination process of the shadowed sea quarks is negligible, the momenta carried by 
each parton component would be conserved separately and approximately, 



< F° /N {Q 2 ) > 2 -< F° c/A {Q 2 ) >2 

< F? /N m >2=< F? /A m > 2 

and the following inequalities will be satisfied for Q 2 > Qq[15] 



c/N(A)\^ J ^l^^ ^ c/N(A) 

As a result, following Re. [O we acquire: 



(27) 



< K$ ( AQ 2 ) >2»< KSiAQ 2 ) >3 (28) 



R s ea(Q , A eff ) = 1 - k s (Q )(Aj ff - 1) = hm « < ^ (gg) ^ (29) 



eff) 



kg(Ql)(A e ff 



H>F*(x,Q§) 



C/n(Qo) >3 



(0? 



>? 



(30) 



where i^ s (q) is the shadowing strength and and depends on the starting scale Ql for the 
evolution to set in. It is easy to verify that the primitive shadowing for gluon in a constituent 
quark is weaker than that for the sea partons[12] at such scale. The form of shadowing 
correction is in terms of evolution kernels is given as: 

p (n 2 a \ _ p i pp. a \ < i71 c/iv(Qg)>3^(Q 2 )+<^/ j vWg)>3^| g (Q 2 )+<^ /i v(Qg)>3^ s (Q 2 ) 

fl>W ' eff) - fi »W0. <F= / ^(Qg)>3^| q (Q 2 )+A<F» / ^(Qg)>3^ 9 (Q2)+i? s (Qg)<F- /JV (Q^)> 3 ^ s (Q2 ) 

(31) 

where K™ b are the evolution kernels. Defining s = InL, these kernel functions are [15]: 



K n qq {Q 2 ) = a n L~< + (1 - a n )L~ a - 



K n NS {Q 2 



a n L~ an + (1 - a n )L~ a ™ - L~ a Ns 



„ _ 7<j<? „n _ 7 :::: " 



(32) 



7 -7 n 

7qg 

7"-7" 1 



III. CONCLUSION 

We have calculated the nucleus structure function ratio in the context of constituent 
quark picture, utilizing the essence of the so called valon model. As we can see from the 
figures; the model describes the data rather well for a wide range of A. It appears that 
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the Q 2 dependence of the ratio of structure functions is week. The shadowing effect is 
quite large for small value of x, (x(0.01) and dies away very rapidly as x increases. Also, 
the shadowing correction becomes more pronounced with increasing atomic number A. We 
further note that in this model, which is based on the constituent quark structure, we 
have not included explicitly the anti-shadowing effects notwithstanding that there are some 
debates on the subject particularly pertinent to the middle range of x. we simply did not 
need it to describe the data. Our anticipation is that such an effect would not be very large 



IV. APPENDIX 

In this appendix we provide numerical values for the coefficients given in equation ( ). 
These are the results of our Q 2 dependence parameterization of parton distributions in a 
constituent quark of proton. They are calculated in the leading order. 

a 8l = -0.119 + 0.0296 exp(s/0.66) 

b sl = 5.510 + 8.50 exp(s/0.95) 
a s2 = -0.080 + 0.066 exp(s/3.12) 
b s2 = -1.194 + 2.303 exp(s/2. 76) 
a gl = -2.490 + 1.027 exp(s/0.69) 

b gl = 5.920 + 6.440 exp(s/0.80) (33) 

a g2 = .174 + .345s - 09 
b g2 = -1.407+ 1.763 exp[s/2.42]; 
a v = -0.300 + 1.254s - 0.368s 2 + 0.034s 3 
b v = 0.196 + 1.286 exp[-s/3. 14] 
C = -0.844 + .451s- 17 
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FIGURES 



FIG. 1. F% function of x at several Q . 



F Al (x Q 2 ) o 

FIG. 2. The ratio R = ' — - for a variety of nuclei at different Q . Dashed-dotted line 
represents the results without the shadowing effects. Others include the shadowing corrections. 
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